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Abstract  -  A  new  Taylor-type  poly  crystalline  model  has  been  developed  to  simulate  the 
evolution  of  crystallographic  texture  and  the  anisotropic  stress-strain  response  during  large 
deformation  of  high  purity  a-titanium  at  room  temperature.  Crystallographic  slip,  deformation 
twinning,  and  slip  inside  the  twinned  regions  were  all  considered  as  contributing  mechanisms  for 
the  plastic  strain  in  the  model.  This  was  accomplished  by  treating  the  dominant  twin  systems  in  a 
given  crystal  as  independent  grains  once  the  total  twin  volume  fraction  in  that  crystal  reached  a 
predetermined  saturation  value.  The  newly  formed  grains  were  allowed  to  independently 
undergo  further  slip  and  the  concomitant  lattice  rotation,  but  further  twinning  was  prohibited. 
New  descriptions  have  been  proposed  for  slip  and  twin  hardening  and  the  complex  coupling 
between  them.  Good  predictions  were  obtained  for  the  overall  anisotropic  stress-strain  response 
and  the  texture  evolution  in  three  different  monotonic  deformation  paths  on  annealed,  initially 
textured  samples  of  high  purity  a-titanium. 
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1.  Introduction 


A  rigorous  formulation  of  the  elastic-plastic  constitutive  relations  for  crystalline  materials 
undergoing  finite  plastic  strains  by  crystallographic  slip  alone  is  now  well  established  in 
literature  [1-8].  However,  many  materials  exhibit  deformation  twinning  as  an  additional  mode  of 
plastic  deformation,  especially  at  low  homologous  temperatures  and/or  high  strain  rates  [9-11]. 
Several  studies  have  reported  extensive  deformation  twinning  in  room  (and  lower)  temperature 
deformation  of  a-Ti  [12-19].  Consequently,  it  is  highly  desirable  to  extend  the  current  crystal 
plasticity  theories  to  include  deformation  twinning  as  an  additional  mode  of  plastic  deformation. 
The  main  obstacle  in  accomplishing  this  goal  is  the  lack  of  an  efficient  method  to  handle  the 
extremely  large  number  of  new  orientations  created  by  deformation  twinning.  Three  different 
approaches  have  been  proposed  in  the  literature  to  address  this  problem: 

(1.)  Predominant  Twin  Reorientation  (PTR)  Method:  This  method  was  originally  proposed  by 
Van  Houtte  [20]  and  improvements  were  made  by  Tome  et  al.  [21].  Staroselsky  and  Anand  [22, 
23]  have  formulated  a  rigorous  and  efficient  numerical  approach  for  the  use  of  rate-independent 
crystal  plasticity  theories  and  applied  it  with  the  PTR  model  to  study  stress-strain  responses  and 
texture  evolution  in  low  stacking  fault  energy  cubic  and  hexagonal  metals.  Kaschner  et  al.  [24] 
employed  the  PTR  method  in  a  visco-plastic  self-consistent  crystal  plasticity  modeling 
framework  to  study  the  role  of  twinning  in  the  hardening  response  of  zirconium  during 
temperature  reloads.  In  the  PTR  scheme,  twinning  is  essentially  treated  as  a  pseudo-slip 
mechanism  while  the  evolution  of  the  volume  fractions  of  the  twinned  regions  is  carefully 
tracked  in  each  grain.  Consequently,  the  twinned  regions  are  not  reoriented  at  the  end  of  each 
time  step.  Instead,  a  statistical  criterion  is  devised  and  employed  to  ensure  that  an  appropriate 
number  of  heavily  twinned  grains  are  completely  reoriented  into  their  predominant  twin 
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orientations;  this  number  is  selected  based  on  the  total  volume  fraction  of  twinned  regions  in  the 
entire  poly  crystal.  A  major  advantage  of  this  method  is  that  the  number  of  total  grain  orientations 
remains  constant  during  the  entire  simulation.  The  disadvantage  of  this  method  is  that  it  can  not 
be  applied  at  the  single  crystal  level  and  is  not  particularly  amenable  to  the  highly  efficient,  total 
Lagrangian,  fully  implicit  time  integration  procedure  developed  recently  for  the  crystal  plasticity 
constitutive  framework  [2], 

(2.)  Volume  Fraction  Transfer  (VFT)  Scheme:  This  approach  was  proposed  by  Tome  et  al.  [21] 
and  employs  weighted  (and  binned)  grain  orientations  to  address  the  problem  of  tracking  the 
large  number  of  new  orientations  created  by  deformation  twinning.  In  this  scheme,  the  relevant 
Euler  space  of  distinct  grain  orientations  (also  called  the  fundamental  zone  of  orientations  [25]) 
is  suitably  binned  and  the  texture  in  the  sample  is  represented  by  the  weights  associated  with  the 
binned  orientations.  Therefore,  only  the  weights  of  the  grain  orientations  need  to  be  suitably 
modified  to  reflect  the  orientation  changes  caused  by  twinning,  and  hence  there  is  no  need  to 
create  new  orientations.  However,  the  disadvantages  of  this  method  are  that  a  very  large  number 
of  bins  are  needed  (especially  for  the  lower  symmetry  hexagonal  metals)  and  that  this  approach 
is  also  not  particularly  amenable  for  implementation  of  the  efficient  fully  implicit  time 
integration  procedure  developed  recently  for  the  crystal  plasticity  constitutive  framework. 

(3.)  Total  Lagrangian  Approach:  Building  on  a  total  Lagrangian  crystal  plasticity  framework  that 
was  initially  formulated  for  materials  that  exhibit  only  crystallographic  slip  [2],  Kalidindi  [26] 
proposed  a  new  interpretation  of  the  multiplicative  decomposition  of  the  total  deformation 
gradient  into  its  elastic  and  plastic  components  when  deformation  twinning  is  added  as  an 
additional  mode  of  plastic  deformation.  The  main  advantages  of  this  method  are  that  it  allows  the 
application  of  the  crystal  plasticity  theory  with  deformation  twinning  to  a  single  crystal  while 
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taking  full  advantage  of  the  efficient  fully  implicit  time  integration  schemes  that  have  been 
previously  developed  and  validated.  The  disadvantage  of  this  method  continues  to  be  the  fact 
that  the  numerical  implementation  of  the  scheme  is  quite  cumbersome  when  slip  inside  the 
twinned  regions  is  to  be  allowed  as  a  significant  contribution  to  the  plastic  deformation  in  the 
crystal.  The  importance  of  including  slip  inside  twins  as  an  additional  mode  of  plastic 
deformation  was  established  in  a  recent  experimental  study  [12,  17,  19],  where  it  was  observed 
that  the  averaged  Taylor  factor  for  the  twinned  regions  in  a  sample  of  high  purity  a-Ti  deformed 
in  simple  compression  was  about  40%  lower  than  that  for  the  matrix  regions,  implying  that  the 
slip  inside  the  twinned  regions  would  be  a  significant  component  of  the  overall  plastic 
deformation  in  this  sample.  In  our  prior  modeling  effort  in  low  stacking  fault  energy  fee  metals, 
the  fact  that  the  potential  twin  systems  are  co-planar  with  the  potential  slip  systems  made  it 
relatively  easy  to  incorporate  slip  inside  twinning  as  a  contribution  to  the  overall  plastic 
deformation  [27],  However,  the  potential  twin  systems  in  hep  metals  are  non-coplanar  with  all  of 
the  potential  slip  systems,  and  therefore  tracking  the  slip  activity  inside  the  already  formed  twins 
becomes  quite  cumbersome  in  this  approach. 

It  is  therefore  clear  that  there  exist  a  number  of  outstanding  problems  in  the  incorporation  of 
deformation  twinning  in  the  current  crystal  plasticity  modeling  framework.  Recent  experimental 
studies  [12,  19]  have  provided  several  new  insights  that,  in  turn,  suggest  new  simplifying 
approximations  to  overcome  some  of  the  difficulties  described  above  in  the  modeling  efforts.  Of 
particular  importance  to  the  work  presented  here  is  the  observation  reported  by  Salem  et  al.  [19] 
that  much  of  the  deformation  twinning  in  high  purity  a-Ti  deformed  in  simple  compression  at 
room  temperature  occurred  in  a  relatively  narrow  moderate  strain  range  of  0.05  to  0.3  (with  most 
of  the  profuse  twinning  actually  occurring  in  an  even  narrower  strain  range  of  0.15-0.20).  At 


4 


higher  strain  levels,  the  grain  structure  fragmented  substantially  and  the  rate  of  deformation 
twinning  decreased  dramatically.  Although  the  precise  reason  for  the  shutdown  of  deformation 
twin  production  at  high  strains  is  not  yet  fully  understood,  it  is  generally  attributed  to  the  much 
higher  stresses  required  to  produce  the  shorter  deformation  twins  in  the  smaller  fragmented 
grains  [28]. 

The  main  objective  of  the  present  paper  is  to  report  on  our  latest  effort  to  develop  an 
improved  Taylor-type  crystal  plasticity  model  for  a-Ti  hep  polycrystals.  This  new  model 
employs  a  grain  fragmentation  concept.  In  this  new  model,  following  the  earlier  approaches  [20, 
29],  deformation  twinning  is  initially  treated  as  a  pseudo-slip  mechanism.  When  a  certain 
volume  fraction  of  any  grain  has  undergone  deformation  twinning,  that  particular  grain  is 
fragmented  into  its  dominant  twin  systems.  In  subsequent  deformation,  neither  the  parent  grain 
nor  any  of  its  newly  formed  fragmented  components  are  allowed  to  produce  deformation  twins. 
The  main  motivation  for  exploring  this  approach  stems  from  the  previously  described 
observations  in  a-Ti  where  it  was  noted  that  deformation  twinning  dominates  the  texture 
evolution  in  the  polycrystal  in  a  short  intermediate  strain  regime  [19,  28].  Therefore,  the 
underlying  approximations  in  the  proposed  approach  should  have  very  little  effect  on  the  texture 
evolution  of  the  polycrystal  at  large  strains.  Furthermore,  since  the  grain  is  allowed  to  fragment 
only  once  and  only  into  its  dominant  twin  components,  the  number  of  new  grains  produced  is 
relatively  small  and  can  be  efficiently  handled  with  only  a  moderate  increase  in  the 
computational  effort.  Finally,  the  proposed  approach  can  be  easily  implemented  in  the  previously 
established,  numerically  efficient,  total  Lagrangian  framework  with  the  associated  fully  implicit 
time  integration  procedure  [2],  In  many  ways,  the  proposed  approach  here  combines  the  best 
aspects  of  the  different  modeling  approaches  described  earlier  [20,  21,  26],  Furthermore,  the  slip- 
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twin  hardening  functions  have  been  expanded  from  the  previously  used  versions  [29],  and  have 
been  observed  to  provide  better  predictions  of  the  anisotropic  stress-strain  responses  in  a  range  of 
deformation  experiments. 

2.  A  New  Crystal  Plasticity  Model  for  High  Purity  a-Ti 

A  new  Taylor-type  crystal  plasticity  model  is  presented  here  for  room  temperature 
deformation  of  high  purity  a-Ti  using  a  notation  that  is  now  standard  in  this  field  [30],  In  this 
notation,  F  represents  the  deformation  gradient  tensor,  L  represents  the  velocity  gradient  tensor, 
and  T  represents  the  Cauchy  stress  tensor.  An  important  feature  of  the  model  presented  here  is 
the  introduction  of  a  grain1  fragmentation  event  based  on  the  accumulated  deformation  twinning 
activity  in  the  crystal.  The  details  of  the  grain  fragmentation  event  and  the  behavior  of  the  grain 
prior  and  subsequent  to  this  fragmentation  event  are  described  next. 

2. 1  Grains  Prior  to  Fragmentation 

The  grains  prior  to  fragmentation  are  modeled  using  essentially  the  approach  described 
previously  by  Kalidindi  [26],  with  some  modifications  in  the  strain  hardening  descriptions.  The 
imposed  total  deformation  gradient  tensor  on  the  crystal  is  multiplicatively  decomposed  into  two 
components  (see  Figure  1): 

F  =  F*FP  .  (1) 


1  The  individual  crystals  comprising  a  polycrystal  are  commonly  referred  to  as  grains  in  metallurgy. 
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Figure  1.  A  schematic  of  the  interpretation  of  the  multiplicative  decomposition  of  deformation 
gradient  when  twinning  is  included  [26]. 


In  this  model,  Fp  denotes  the  plastic  deformation  gradient  tensor  that  describes  the  overall 
(effective)  shape  change  induced  in  the  grain  as  a  consequence  of  the  crystallographic  slip  and 
deformation  twinning  in  the  grain,  while  F*  denotes  the  additional  transformation  (including 
elastic  stretch  and  rotations)  needed  to  account  for  the  total  imposed  shape  change  induced  by  F. 
It  is  further  assumed  that  the  slip  processes  leave  the  lattice  orientation  unchanged,  while 
deformation  twinning  rotates  the  lattice  into  a  priori  defined  orientation  (for  the  hexagonal 
crystals  of  interest  in  this  paper  this  relationship  is  conveniently  described  as  a  180  degree 
rotation  about  the  twin  plane  normal  [9]).  Consequently,  the  orientations  of  the  twinned  and 
untwined  regions  of  the  grain  in  the  hypothetical  intermediate  configuration  shown  in  Figure  1 
(after  the  application  of  Fp )  are  fully  known  a  priori,  given  the  initial  orientation  of  the  grain.  As 
a  consequence  of  the  above  definition  ofFp,  the  polar  decomposition  of  F*  provides  the  elastic 
stretch  and  the  overall  lattice  rotation  between  the  intermediate  configuration  and  the  deformed 
configuration.  In  Figure  1,  only  one  twin  system  is  shown  for  clarity.  However,  multiple  twin 
systems  are  allowed  in  one  grain.  Furthermore,  the  twinned  region  belonging  to  one  twin  system 
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is  idealized  as  a  continuous  block,  while  in  reality,  deformation  twins  occur  as  thin  plates.  This 
simplification  was  motivated  by  the  reported  experimental  observations  in  literature  [12,  17,  31], 
which  suggest  that  the  deformation  twins  of  a  particular  twin  system  occur  as  clusters  of  roughly 
parallel  plates  of  similar  lattice  orientation  in  any  given  grain.  The  fact  that  the  twins  occur  as 
clusters  of  parallel  plates  does  have  a  major  influence  on  the  hardening  response  of  the  crystal, 
and  this  aspect  will  be  treated  later  using  appropriate  phenomenological  slip  and  twin  hardening 
laws. 

The  constitutive  equation  for  the  elastic  response  of  a  single  crystal  is  expressed  as 

T*=c[e‘],  (2) 

where  C  represents  the  fourth-order  anisotropic  elasticity  tensor,  T*  and  E*  are  a  pair  of  work 
conjugate  stress  and  strain  measures  defined  as 

T*  =F*~1{(detF*)T}F"T'  ,  (3) 

E*  =I{F*rF*  -1}.  (4) 

The  evolution  of  plastic  deformation  gradient  can  be  expressed  as 

Fp  =  LPFP ,  (5) 

where  Lp  is  plastic  velocity  gradient  tensor  expressed  as 

Ns  Nm 

L  p=ZrS:+S/"^Sf.  (6) 

a  p 

The  two  terms  on  the  right  hand  side  of  Eq.  (6)  represent  the  contributions  to  plastic  deformation 
by  slip  and  deformation  twinning  respectively.  S  denotes  the  unit  slip  (twin)  tensor,  defined  as 
the  dyadic  product  of  two  orthogonal  unit  vectors  denoting  the  slip  (twin  shear)  direction  and  the 
slip  (twin)  plane  normal,  respectively.  The  subscript  o  on  S  reminds  us  that  these  are  defined 
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using  the  initial  crystal  orientation  and  therefore  the  S  tensors  are  known  a  priori.  The 
superscripts  a  and  p  enumerate  the  available  slip  systems  (total  of  Ns )  and  the  available  twin 
systems  (total  of  M1"' ),  respectively.  ya  represents  the  slip  shear  rate  on  slip  system  a.  Unlike 

slip,  deformation  twinning  is  characterized  by  a  constant  amount  of  shear  and  the  volume 
fraction  of  the  crystal  experiencing  deformation  twinning  evolves  with  the  imposed  deformation. 
Treating  deformation  twinning  as  a  pseudo-slip  mechanism,  the  homogenized  shear  rate  in  the 

crystal  is  described  by  ,  where  yfw  represents  the  amount  of  constant  shear  associated  with 

twin  system  p  and  /^denotes  its  volume  fraction  in  the  given  crystal. 


Power-law  relations  have  been  employed  to  quantify  the  plastic  shearing  rate  on  slip  and  twin 
systems  using  the  visco-plastic  approach  of  Asaro  and  Needleman  [4]: 


ra  =Y0 


a 


a 


m 


sign(Ta), 


(7) 


Yo 

Tp 

/ 

/  tw 

s13 

r 

o 

if 

if 


:fl  >0 


<0 


(8) 


In  this  study,  the  rate  sensitivity  parameter,  m,  was  assumed  to  be  the  same  for  both  slip  and 
twinning.  A  very  low  value  of  m  =  0.02  was  used  to  simulate  almost  rate-independent  behavior 
of  metals  at  room  temperature.  The  reference  slip  rat e,y0,  was  arbitrarily  set  as  0.001  s'1  to 


reflect  our  interest  in  quasi-static  loading  conditions.  ra  and  s a  represent  the  resolved  shear 
stress  and  the  shear  resistance  for  a  particular  slip  or  twin  system.  Note  that,  unlike  slip,  a 
positive  resolved  shear  stress  is  needed  for  twinning.  The  resolved  shear  stress  for  both  slip  and 
twin  systems  can  be  defined  as  [2] 
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Description  of  the  evolution  of  the  slip  and  twin  resistances  during  plastic  deformation  has 
been  a  very  difficult  problem  in  the  development  of  robust  crystal  plasticity  models  for  hep 
metals  such  as  Ti.  The  slip-twin  interactions  are  fairly  complex  [12,  17,  19,  32,  33],  and  there  is 
thus  far  only  a  limited  amount  of  quantitative  experimental  data  available.  In  this  study, 
following  the  current  efforts  in  literature,  we  employed  a  phenomenological  description  of  the 
slip  and  twin  hardening  laws.  Some  of  the  most  successful  phenomenological  descriptions  to 
date  have  been  the  saturation-type  hardening  laws  that  can  be  generically  expressed  as 
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where  h"  and  s“  represent  the  hardening  rate  and  the  saturation  value  associated  with  the  slip 


system  a,  respectively.  Implicit  in  Eq.  (10)  is  the  assumption  of  isotropic  latent  hardening  which 
implies  that  activity  on  one  slip  system  results  in  equal  hardening  of  all  potential  slip  systems.  In 
a  recent  paper,  Salem  et  al.  [29]  used  extended  versions  of  the  saturation-type  hardening 
functions  to  capture  the  complex  interactions  between  slip  and  twinning.  In  that  study,  Eq.  (10) 
has  been  used  for  slip  hardening,  while  allowing  the  slip  hardening  rate  and  the  saturation  values 
to  evolve  with  deformation  twin  activity  as 

K  =h,(\  +  c(£fy),  (11) 
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The  functional  forms  of  Eqs.  (11)  and  (12)  were  motivated  by  experimental  observations  in 
deformation  studies  on  high-purity  a-Ti  [12,  17,  19,  28].  The  parameters  hs  and  Sso  denote  the 
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hardening  rate  and  the  saturation  value,  respectively,  in  the  absence  of  twinning.  C,  s  and  b 
are  hardening  parameters  that  aim  to  capture  phenomenologically  the  complex  interactions 

between  slip  and  deformation  twinning.  ^  fp  denotes  the  total  twin  volume  fraction  in  the  grain. 

P 


In  the  present  work,  we  found  it  essential  to  allow  the  hardening  parameters  in  Eqs.  (11)  and 
(12)  to  take  on  different  values  for  the  different  slip  families  in  order  to  obtain  better  predictions 
of  the  anisotropic  stress-strain  responses  in  a-Ti  subjected  to  different  deformation  paths.  Table 
1  provides  a  summary  of  the  different  slip  and  twin  families  considered  in  this  study.  A  total  of 
18  slip  systems  belonging  to  three  distinct  slip  families  and  12  twin  systems  belonging  to  two 
distinct  twin  families  were  considered.  Consequently,  Eqs.  (10)  -  (12)  were  reformulated  for  the 
present  study  as 
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Tablel.  Summary  of  Slip  and  Twin  systems 


Slip  Systems 

Twin  Systems 

Family 

Prismatic  <a> 

Basal  <a> 

Pyramidal<c+a> 

Compressive 

Tensile 

Plane  &  Direction 

{1  ol 0}  <  1 120  > 

{0001}  <1120: 

{10T 1}  <  1 123  > 

{1122}  <  1123  > 

{l  ol  2}  <  Toi  l  > 
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Number  of  systems 

3 

3 

12 

6 

6 

Initial  resistance 

„Pri 

bas 

A 0 

Qpyr 

*o 

„  com 

Stw 

ten 

Atw 

Slip  family  specific 
hardening  parameters 

ljPri  cPri 

nso  ?  **  SO 

ibas  bas 

nSO  5  ^  SO 

fopyr  pyr 
nso  ?  ^  so 

No  hardening  parameters 

Common  hardening 
parameters 

C,  b,  s ' Spr 

needed  for  twin  systems. 

The  experimental  observations  of  Salem  et  al.  [12,  17]  indicate  that  twinning  occurs  profusely 
in  a  short  intermediate  strain  range.  Consequently,  we  have  decided  to  treat  the  twin  resistances 
as  constants  prior  to  grain  fragmentation. 

2.2  Criterion  for  Grain  Fragmentation 

In  the  present  model,  when  the  twin  volume  fraction  in  a  given  grain  (^/^  )  reaches  a 

p 

predetermined  saturation  value  (  fsat ),  that  grain  is  fragmented  into  parts:  a  parent  grain 
(corresponding  to  the  untwinned  or  matrix  region)  and  several  offspring  grains  (corresponding  to 
the  dominant  twinning  activity  in  that  grain).  In  this  study,  fsat  was  set  equal  to  0.4  based  on  the 

experimental  observations  reported  by  Salem  et  al.  [12,  19].  The  volume  fractions  of  the  parent 
and  the  offspring  grains  are  determined  by  the  accumulated  twinning  activity  in  the  grain  up  to 
the  point  of  fragmentation.  It  is  ensured  that  the  sum  of  the  volume  fractions  of  the  parent  and 
the  offspring  grains  is  equal  to  one.  Once  the  grain  is  fragmented,  the  volume  fractions  of  the 
newly  created  grains  are  kept  constant  during  the  rest  of  the  deformation.  Consequently,  the 
number  of  new  orientations  produced  in  this  model  simply  equals  to  the  number  of  offspring 
grains  created  in  the  fragmentation  process.  In  this  study,  only  the  twin  systems  with  a  volume 
fraction  of  over  0.1  were  represented  in  the  newly  created  grains.  The  volume  fractions 
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associated  with  the  non-dominant  twin  systems  (defined  here  as  fp  <  0.1 )  have  been  transferred 
to  the  most  dominant  twin  system. 

2.3  Grains  after  Fragmentation 

In  this  study,  the  slip  resistances  of  the  various  slip  systems  in  the  parent  and  offspring  grains 
have  been  assigned  the  same  values  that  were  in  the  grain  prior  to  fragmentation.  However,  the 
resistance  for  any  twin  system  in  the  parent  and  the  offspring  grains  is  assumed  to  be  very  high 
so  that  no  further  twinning  can  be  activated.  The  evolution  of  slip  resistances  in  subsequent 
deformation  continues  as  described  by  Eqs.  (13)  -  (15).  Note  that  the  slip  activities  in  the 
offspring  grains  are  tantamount  to  slip  in  the  twinned  regions  of  the  original  grain.  After  grain 
fragmentation,  the  newly  formed  grains  (parent  and  offspring  grains)  are  essentially  treated  as 
independent  grains.  Moreover,  further  plastic  deformation  in  these  grains  is  assumed  to  be  fully 
accommodated  by  crystallographic  slip  alone.  Note  that,  in  the  total  Lagrangian  framework 
employed  in  this  work,  the  slip  tensors  (S“in  Eq.  (6))  in  the  newly  created  grains  are  defined 

based  on  the  orientations  of  the  twinned  regions  in  the  intermediate  relaxed  configuration.  As 
noted  earlier,  the  twin  orientations  in  the  intermediate  configuration  are  known  a  priori,  and  are 
conveniently  described  by  a  180  degree  rotation  about  the  twin  habit  plane  for  the  a-Ti  single 
crystals  studied  here. 

The  treatment  of  the  grains  before  and  after  fragmentation  is  summarized  in  Figure  2.  The 
following  salient  features  of  the  new  model  presented  here  are  worth  noting: 

(1.)  Before  fragmentation,  deformation  twinning  is  treated  as  a  pseudo-slip  mechanism.  So  there 
is  only  a  single  decomposition  of  the  total  deformation  gradient  tensor.  Consequently,  the 
matrix  and  the  twinned  regions  have  a  pre-defmed  orientation  relationship.  In  other  words, 
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the  twin  regions  rotate  with  the  matrix  as  if  they  were  rigid  inclusions  attached  to  the  matrix 
with  a  specific  orientation  relationship. 

(2.)  After  fragmentation,  the  newly  formed  grains  behave  independent  of  the  parent  grain  and 
each  other.  Following  the  Taylor  approach  [34],  the  matrix  and  twin  are  assumed  to 
experience  the  same  total  deformation  gradient  tensor.  However,  the  plastic  deformation 
gradient,  the  elastic  deformation  gradient,  and  the  lattice  rotation  are  expected  to  differ 
substantially  among  the  parent  and  each  of  the  offspring  grains.  Note  also  that  the  special 
orientation  relationship  between  the  parent  and  offspring  grains  is  still  maintained  in  the 
intermediate  relaxed  configuration  (even  after  grain  fragmentation),  but  is  expected  not  to  be 
retained  in  the  final  configuration  (because  the  parent  and  offspring  grains  are  generally 
expected  to  experience  different  lattice  rotations  after  the  fragmentation  event). 
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Figure  2.  Summary  of  the  multiplicative  decomposition  of  the  total  deformation  gradient 
implemented  in  the  model  presented  in  this  study,  pg  denotes  the  parent  grain  and  og  denotes  the 
offspring  grain.  Only  one  offspring  is  shown  in  the  figure  for  clarity. 

The  response  of  a  poly  crystalline  aggregate  is  obtained  here  using  the  extended  Taylor 
assumption  of  uniform  deformation  gradient  tensor  in  the  entire  polycrystal  [34],  Consequently, 
the  overall  Cauchy  stress  tensor  in  the  polycrystal  is  assumed  to  be  given  by  the  volume 
averaged  value  of  the  Cauchy  stress  tensor  in  all  of  the  constituent  grains  (including  all  of  the 
parent  and  the  offspring  grains). 

3.  Experiments  on  High  Purity  a-Ti 

The  material  used  in  this  study  was  a  high  purity  a-titanium  (99.9998%)  supplied  by  the  Alta 
Group  of  Johnson  Matthey  Electronics,  Inc.,  Spokane,  WA.  The  received  material  was  a  circular 
plate  with  352  mm  diameter  and  12  mm  thickness.  This  material  was  recrystallized  at  800°C  for 
1  hour  and  followed  by  water  quench,  producing  an  equiaxed  grain  structure  with  a  30  pm 
average  grain  size.  The  annealed  sample  showed  a  strong  fiber  texture  (see  Figure  5).  The  (0001) 
pole  figures  in  figure  5  indicate  that  the  c-axes  of  many  grains  in  this  sample  lie  about  20-35 
degrees  to  the  plate  normal  (labeled  ND).  Directions  RD  and  TD  were  marked  arbitrarily  on  the 
plane  of  the  plate  surface  as  shown  in  Figure  5.  Additional  information  about  this  material  can  be 
obtained  from  the  references  [12,  17,  19]. 

Three  stress-strain  measurements  (two  from  simple  compression  tests  and  one  from  simple  shear 
test)  on  this  material  and  the  corresponding  strain  hardening  plots  were  available  from  previous 
studies  [12,  17,  29].  Measurements  of  the  initial  texture  (Figure  5)  and  the  measurement  of  deformed 
texture  in  simple  shear  test  were  also  available  from  the  earlier  work  [29],  Two  new  measurements 
of  the  deformed  textures  in  simple  compression  tests  were  acquired  for  this  study.  These  new 
measurements  were  obtained  by  X-ray  reflection  technique  on  a  Scintag  XI  5-axis  pole  figure 
goniometer  using  Cu-Ka  radiation.  The  beam  is  point  source,  circularly  collimated  to  0.8mm.  The 
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detector  is  a  Peltier  cooled  solid  state  detector  with  a  slit  of  2mm.  The  samples  were  a  minimum  of  1 
cm2.  The  (10 1 0)  ,(0001)  ,  (1011)  ,  (1012)  ,  (1120)  and  (10 13)  pole  figures  were  measured.  The 
samples  were  oscillated  4  mm  to  collect  data  from  as  many  grains  as  possible.  Data  was  collected 
out  to  80  sample  tilt  in  5  sample  rotation  and  tilt  increments  with  a  two  second  counting  time.  The 
data  was  analyzed  using  the  popLA  software  [35],  to  produce  corrected  and  recalculated  pole  figures 
and  orientation  distributions.  The  raw  data  was  corrected  for  background  and  defocusing  using  a 
correction  file  generated  from  a  titanium  sample  with  a  random  texture.  The  data  was  then  run 
through  a  harmonics  algorithm  to  extrapolate  the  outer  fringes  of  the  pole  figure  data  and  then  the 
data  was  re-normalized.  This  data  was  run  through  the  WIMV  algorithm  to  produce  recalculated 
pole  figures  and  orientation  distributions. 

4.  Calibration  and  Evaluation  of  Model 


4.1  Calibration:  determination  of  model  parameters 

There  are  a  total  of  fourteen  material  parameters  in  the  crystal  plasticity  model  presented  in 
the  previous  section  (see  Table  1).  Our  goal  is  to  establish  the  values  of  these  parameters  by 
curve-fitting  the  predicted  stress-strain  responses  in  selected  deformation  modes  to  the 
corresponding  experimental  measurements.  In  previous  studies  [29],  we  have  established  an 
approach  for  accomplishing  this  task.  Our  approach  recognizes  the  influence  the  different 
parameters  have  on  the  form  of  the  predicted  stress-strain  curves  and  uses  repeated  trials  until  the 
predictions  match  the  measurements.  In  this  process,  in  the  present  study,  we  noted  that  some  of 
the  parameters  for  the  different  slip  families  in  Eqs.  (13)-(15)  exhibited  values  very  close  to  each 
other.  Based  on  these  observations,  we  made  the  following  assumptions  for  the  present  study 
that  reduced  the  hardening  parameters  from  fourteen  to  eleven: 


^  pi'i  —}ibas  —  j^Pri~bas 


fas  _  pyr  _  bas-pyr 
** so  ^  so  ^  so  ? 


com  _  ten  _ 

^tw  ^ tw  ^ tw 


(16) 


The  values  of  all  the  eleven  parameters  were  established  by  fitting  the  predicted  stress-strain 
curves  to  the  corresponding  measurements  in  two  different  monotonic  deformation  modes 
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(simple  compression  and  simple  shear).  Briefly,  our  strategy  for  establishing  the  hardening 
parameters  comprised  of  the  following  steps  [29]: 

pri  has  py y 

(1)  The  values  of  initial  resistances  ( S0  ,  so  and  S0  )  on  different  slip  families  (  basal  <a>, 

prism  <a>  and  pyramidal  <c+a>)  were  determined  by  fitting  the  predicted  yield  strengths  in 
simple  compression  along  ND  and  in  simple  shear  to  the  corresponding  measurements  (see 
Figure  3). 

(2)  The  value  of  initial  resistance  ( stw )  for  twinning  was  determined  by  matching  the  start  of 

stage  B  in  the  predicted  strain  hardening  plot  in  simple  compression  along  ND  (Point  2  in  Figure 
4)  to  the  corresponding  measurement.  In  a  recent  report  [12],  this  point  on  the  strain  hardening 
plot  curve  has  been  correlated  with  the  onset  of  deformation  twinning. 

(3)  The  values  of  slip  hardening  parameters  ( hpr0l~bas ,  hpyr ,  sp''  and  sbsa0s~pyr )  were  determined  by 

fitting  the  stage  A  of  strain  hardening  plot  (see  Figure  4)  and  the  stress-strain  curve  in  simple 
shear  (see  Figure  3)  where  the  volume  fraction  of  deformation  twinning  is  not  significant  and 
therefore  slip  is  considered  as  the  main  mechanism  to  accommodate  plastic  deformation. 

(4)  The  values  of  twin  hardening  parameters  (b  and  C)  were  determined  by  fitting  the  stage  B  of 
strain  hardening  plot  (see  Figure  4)  where  the  increment  of  strain  hardening  rate  is  correlated  to 
the  increasing  twin  volume  fraction. 

(5)  The  value  of  twin  hardening  parameter  ( s  )  was  determined  by  fitting  the  saturated  flow 
stress  in  simple  compression  along  ND  (see  Figure  3). 
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Figure  3.  Comparison  of  predicted  (P)  and  measured  (M)  equivalent  stress-equivalent  strain 
curves  for  different  deformation  modes  on  high  purity  a-Ti. 


Figure  4.  Strain  hardening  response  of  a-Ti  in  simple  compression  along  ND.  The  ordinate  is  the 
normalized  slope  of  the  stress-strain  curve  and  the  abscissa  is  the  normalized  plastic  flow  stress, 
and  G  is  theoretical  shear  modulus. 

Using  above  strategy,  the  values  of  all  model  parameters  were  established  (see  Table2). 
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Table2.  Summary  of  Parameter  Values 


Model 

parameters 

?Pri 

bas 

% 

„pyr 

Stw 

7  pri-bas 
nso 

hPyr 

nso 

Values 

30MPa 

150MPa 

120MPa 

125MPa 

15  MPa 

300MPa 

Model 

parameters 

?Pri 

^  SO 

bas- pyr 

Sso 

Spr 

C 

b 

Values 

lOOMPa 

300MPa 

lOOMPa 

25 

2 

3.2  Evaluation:  comparison  with  measurements 

In  this  section,  the  proposed  model  and  hardening  laws  were  evaluated  quantitatively  by 
directly  comparing  the  predictions  to  the  measurements  that  were  not  used  in  the  calibration 
process.  The  predicted  stress-strain  response  in  simple  compression  along  TD  showed  excellent 
agreement  with  the  measurement  (see  Figure  3).  It  is  worth  noting  that  the  starting  texture  for  the 
compression  test  along  TD  is  significantly  different  compared  to  the  starting  texture  for 
compression  test  along  ND  (see  Figure  5).  Note  also  that  the  flow  stress  in  compression  along 
TD  is  substantially  lower  than  the  flow  stress  in  compression  along  ND  (see  Figure  3).  The 
predicted  stress-strain  response  for  simple  shear  also  showed  good  agreement  with  measurement 
up  to  strain  y  «  0.6 ,  but  an  overestimation  of  about  20%  on  flow  stress  is  observed  at  large  strain 
( y  ~  1  -0 ;  see  Figure  3).  It  is  worth  noting  that  at  large  shear  strains  the  shear  samples  often 
develop  macroscale  shear  bands  and  this  could  lower  the  measured  flow  stresses. 

The  stress-strain  curves  in  Figure  3  were  plotted  as  equivalent  stress  -  equivalent  strain  curves. 
While  the  concept  of  a  Mises-equivalent  stress  has  no  relevance  for  anisotropic  plasticity,  it 
helps  us  understand  the  degree  of  anisotropy  exhibited  by  the  material.  If  the  material  exhibited 
isotropic  plastic  response,  all  of  the  curves  shown  in  Figure  3  should  be  identical.  The  spread 
between  the  stress-strain  responses  in  the  different  modes  (about  a  factor  of  2  between  simple 
compression  and  simple  shear)  in  this  figure  provides  us  a  measure  of  the  degree  of  anisotropy 
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exhibited  by  the  high  purity  a-Ti  that  is  largely  attributable  to  the  underlying  crystallographic 
texture  in  the  sample.  This  degree  of  anisotropy  appears  to  be  well  captured  by  the  Taylor-type 
crystal  plasticity  model  presented  in  this  paper. 


ND 

000  1 


RD 


Figure  5.  Comparison  of  initial  texture  for  compression  tests  on  a-Ti  along  ND  and  TD.  The 
compressed  samples  were  cut  from  a  circular  plate. 

Texture  predictions  for  simple  compression  along  ND  at  two  strains  (  s  =  -0.22  and 
£•  =  -1.00)  showed  excellent  agreement  with  the  measurements  (see  Figure  6  and  Figure  7). 
At£'  =  -0.22,  both  the  prediction  and  the  measurement  exhibit  two  major  c-axis  fiber  texture 
components  (easily  seen  in  (0001)  pole  figure).  The  first  fiber  has  its  c-axis  about  15-30°  from 
ND,  while  the  second  fiber  has  its  c-axis  about  75-90°  from  ND.  The  first  component  was 
associated  with  matrix,  and  the  second  component  with  twinning,  in  a  prior  study  [19].  A  similar 
texture  is  observed  in  (0001)  pole  figure  at£-  =  -1.00,  where  the  intensity  of  the  component 
associated  with  twinning  increased  substantially. 


Measurement 


20 


Prediction 


1010  0001  1011 


Figure  6.  Comparison  of  predicted  and  measured  textures  at  s  =  -0.22  in  simple  compression  of 
a-Ti  along  ND. 
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Figure  7.  Comparison  of  predicted  and  measured  textures  at  s  =  -1.00  in  simple  compression  of 
a-Ti  along  ND. 


Texture  predictions  for  simple  shear  at  y  =  -1.00  showed  reasonable  agreement  with  the 


corresponding  measurement  (see  Figure  8),  especially  in  the  (1010)  pole  figure  where  six 


strong  texture  components  are  clearly  seen.  The  predicted  (0001)  pole  figure  captured  the  strong 
texture  component  located  about  20-40°  from  ND.  However,  the  predicted  (0001)  pole  figures 


missed  the  weaker  texture  components  located  around  the  rim  (see  Figure  8). 
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Figure  8.  Comparison  of  predicted  and  measured  textures  at  y  =  -1.00  in  simple  shear  of  a-Ti. 


3.3  Discussion 

The  poly  crystal  plasticity  model  presented  here  is  built  on  the  Taylor  assumption  of  uniform 
deformation  gradient  in  all  of  the  constituent  grains.  It  should  be  intuitively  expected  that  this 
gross  simplification  should  have  a  strong  effect  on  the  predictions.  In  fact,  it  should  be  expected 
that  the  effect  of  the  Taylor  assumption  would  be  stronger  on  the  hep  metals  studied  here 
compared  to  the  more  plastically  isotropic  cubic  metals  studied  in  previous  work  [1,  2,  8,  1 1,  27]. 
It  is  our  intuition  that  many  of  the  discrepancies  reported  here  between  the  measurements  and  the 
predictions  are  all  attributable  to  the  Taylor  assumption  employed  in  this  study.  We  are  currently 
attempting  to  implement  the  crystal  plasticity  model  described  here  in  a  finite  element 
framework  which  averts  the  need  for  the  simplifying  Taylor  assumption. 

We  also  note  that  only  monotonic  deformation  paths  have  been  studied  thus  far.  There  is 
therefore  a  clear  need  to  extend  the  studies  described  here  to  deformation  path  change 
experiments.  It  is  anticipated  that  such  studies  will  help  us  improve  further  the  hardening 
descriptions  [36-38], 

4.  Conclusions 

A  new  crystal  plasticity  model  has  been  formulated  to  simulate  the  anisotropic  stress-strain 
response  and  texture  evolution  for  a-titanium  during  large  plastic  strains  at  room  temperature. 
The  major  new  features  of  the  model  include:  (a)  incorporation  of  slip  inside  twins  as  an 
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significant  contributor  to  accommodating  the  overall  imposed  plastic  deformation,  and  (b) 
extension  of  slip  and  twin  hardening  laws  to  treat  separately  the  hardening  behavior  of  the 
different  slip  families  (prismatic<a>,  basal<a>,  and  pyramidal  <c+a>)  using  hardening 
parameters  that  are  all  coupled  to  the  extent  of  deformation  twinning  in  the  sample.  Reasonable 
agreement  between  model  predictions  and  the  measurements  has  been  observed  for  both  the 
anisotropic  stress-strain  responses  and  the  deformation  texture  evolution  in  three  different 
monotonic  deformation  paths:  (1)  simple  compression  along  ND,  (2)  simple  compression  along 
TD,  (3)  simple  shear  in  the  RD-TD  plane. 
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